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PACS 89 . 75 . Da - Systems obeying scaling laws 

Abstract. - This Letter studies the critical point as well as the discontinuity of a class of explo- 
sive site percolation in Erdos and Renyi (ER) random network. The class of the percolation is 
implemented by introducing a best-of-m rule. Two major results are found: i). For any specific 
771, the critical percolation point scales with the average degree of the network while its exponent 
associated with m is bounded by —1 and ^ —0.5. ii). Discontinuous percolation could occur 
on sparse networks if and only if m approaches infinite. These results not only generalize some 
conclusions of ordinary percolation but also provide new insights to the network robustness. 



^^ Percolation describes the connectivity of a graph by con- 
^^tmuously occupying links or nodes. Due to its wide appli- 
L_iation in a variety of area such as epidemics, nuclear mul- 
tifragmentation and network robustness etc., it has been 
\^an active subject of research for decades [T]-|6]. Perhaps of 
JiT t'he greatest importance in percolation studies is its order 
^^of phase transition and the location of critical point. In 
fS| almost all cases percolation on graphs is shown to be con- 
^Dtinuous (higher than first order) while critical points for 
^^f^ networks of various structure and dimension have been de- 
^^ ter mined. In an ER random network with average degree 
CN k, as a well studied case, ordinary site percolation causes 
\ ]a continuous phase transition at critical point tc{k) = 1/k 



k>H ■ Recently a new kind of percolation, named explosive 
^ percolation was proposed. By introducing a proper com- 

_C^petitive mechanism, it was first found by Achlioptas, 
D'Souza, and Spencer and was subsequently studied in- 
tensely by other scientists that the bond percolation in 
random networks could be discontinuous J9HT7]. However, 
further numerical and theoretical studies demonstrated 
that such percolation is actually continuous in the ther- 
modynamic limit but has unusually small critical expo- 
nent [T8H2T]. But this is not the end of the story. For 
instance, a very recent study pointed out that the be- 
havior of the explosive percolation transition depends on 
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detailed dynamic rule [22]. When dynamic rules are de- 
signed to suppress the growth of all clusters, the explosive 
percolation transition could be discontinuous. Another re- 
cent study on explosive site percolation on square lattice 
also claimed the existence of discontinuous phase ] transi- 
tion [23]. Therefore whether the explosive percolation is 
indeed discontinuous or continuous is still controversial. 

Despite the extensive studies on discontinuity of explo- 
sive percolation, the other important property, namely the 
location of the critical point, has not been studied sys- 
tematically. Previous studies presented the related results 
only for their special models, which neither provide any 
general conclusions nor help to understand its physical 
meaning in network dynamics. In this Letter, critical per- 
colation point is the prior issue to discuss. This is partly 
motivated by some earlier works on network robustness 
under attack [231[25]. While most of their attack strate- 
gies are based on degree information, real-world networks 
may suffer more diverse structure frangibility to which a 
practical attack strategy could be more specific. The con- 
cept of explosive percolation allows us to explore these 
questions in a general way and could provide new insights 
to network robustness. One of the aim of the Letter is 
to show that in addition to discontinuity or not, explosive 
percolation has other appealing properties worth study- 
ing. The discontinuity of the percolation and the related 
critical exponents are also discussed properly. 

Let us begin with our explosive rule. Consider a con- 
nected ER random network with total number of nodes 
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Fig. 1: Visualization of the explosive rule for the case ofm = 2. The 
empty nodes are unoccupied nodes. Four circles of different color 
and the corresponding inside nodes represent four occupied clusters. 
Two golden nodes (Vi and V2) are the candidates selected randomly 
and compete for occupation. The one connecting the clusters of 
the smallest total size will be chosen. In this case, the node V2 is 
finally chosen since it connects the clusters of total size 3, smaller 
than 5, the total size of the clusters that Vi connects. Note that the 
candidates are selected from all possible empty nodes. If a selected 
candidate is an isolated single node of ER network, this node itself 
represents a cluster of size equaling to 1. 



TV and average degree k. At each time step t, defined as 
the fraction of number of occupied nodes of N^ m empty 
nodes are selected randomly as candidates to be occupied, 
but only the one which minimizes the sum of the size of 
the clusters that itself connects is finally chosen, as visu- 
alized in Fig [TJ If there is more than such a node, we 
choose one of them randomly. This competitive process is 
the so-called best- of -m rule which has also been applied to 
bond percolation [T6l[T9]. Repeating this process, the sys- 
tem will eventually become percolated at a critical point 
tc = tc{m,k^N), while its critical behavior depends cru- 
cially on the value of m. With the increase of m, the 
growth of the order parameter 5, defined as 5 = Si/N 
where ^i is the the size of the largest connected clus- 
ter, is effectively suppressed, leading to gradually delayed 
critical points and a discontinuous-hke jump in its crit- 
ical dynamics (Fig. [2](a)). For the case of m = 1, the 
traditional continuous site percolation is recovered while 
for 771 = A^ the system becomes the most explosive. In 
most cases, percolation concerns the critical behavior in 
the thermodynamics limit, so we focus on critical point 
in the thermodynamics limit tc{m^ /c), which is defined as 
tc{m^ k) = tc{m^ k^N ^ 00). 

Further inspection on Fig. [2](a) reveals that the in- 
crease rate of the critical point decreases with m. In 
fact, the increase rate vanishes so quickly that the crit- 
ical point is conjectured to have a nontrivial limitation 
Tc{k) = tc{m = N ^ 00, k) < 1 . This conjecture can be 
reached by contradiction. Suppose such limitation does 
not exist, i.e. Tc = 1, it is expected that just before this 
Tc, the network has only N^ {r] < 1) unoccupied nodes 
and N^~^{a < 1) isolated occupied clusters of similar size 




Fig. 2: (a): Site percolation process for the case of /c = 6. With 
the increase of m, critical point becomes much more delayed, but its 
delayed rate decreases rapidly. Eg., the delayed time is 0.198 between 
m = 1 and 2 but only 0.038 between m = 10 and 20. The change of 
m is ten times larger but the delayed rate is five times smaller than 
the previous one. The figure is a single simulation for network size of 
5 X 10^ nodes, (b): The increase of tc with m for the case oi k = 6. 
Blue solid point is the simulated data while red dashed line is the fit 
by using Eq.(2). In this case, Tc = 0.633, /(m) = -0.490m-0-'^94 
Inset: the power-law decay of f'(m) obtained by calculating the 
difference oi tc(m) for successive m from 1 to 16. In this simulation, 
numerical calculations for tc(m, k, N) are carried out up to network 
size of 10^ nodes and are averaged over 50 realizations. 



of 0{N^) as global competition in explosive bond perco- 
lation. [19]. Denoting k' the average degree within each 
of the cluster, these clusters surely belong to the ensem- 
ble of subgraph of N^ nodes and k'N'^/2 links. On the 
other hand in ER random network, the average number of 
appearance G of such subgraph can be calculated as 



(G)-7V(^-")^"(^-^^ 



(1) 



Consistently we should have (G) > N^~^, leading to 
k' <2. However k' >2 must hold because each occupied 
cluster is connected. Then in the thermodynamic limit it 
comes to the only possible solution k' ^ 2. While for a 
connected random graph, the average degree k is gener- 
ally larger than 2, so for each of the isolated cluster there 
are about (k — 2)N^ links stem from the cluster to the 
remaining unoccupied nodes, leading to 0{N) such links 
totally. We thus arrive at a contradiction since the unoc- 
cupied nodes can only provide kN^ /2 < 0{N) links. Note 
that although the proof is made for /c > 2, it does not rule 
out the possibility of its validation for /c < 2. Actually 
Tc = 1 only occurs when /c ^ 1 as will be discussed later. 
It is then nontrivial to find such limitation Tc. Defining 
the susceptibility as x = A^\/(5'^) — (5')^, which quan- 
tifies the amplitude of the fiuctuations of the size of 
the largest cluster, the critical point tc{m,k,N) of fi- 
nite network size N is said to locate at which x reaches 
its maximum. Then with a finite size scaling function 
tc{m,k,N) — tc(rn,k) ~ 7V~^/^, the critical point in the 
thermodynamic limit tc{m, k) can be precisely determined 
[12]. By this method, tc{m, k) for each m = 2, 3, 5, 10, 20 
and each /c = 4, 5, 6, 8, 10, 14 is calculated as a preparation 
for finding Tc. Tc for a particular k can be written as 

tc{m)=Tc^f{m). (2) 

It is valid that tc(m+l)—tc(m) = /(m+1) — /(m) ^ f'{m) 
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, which indicates that the formula of f{m) can be de- 
rived by calculating the difference of tdrn). Then fur- 
ther calculation of tc{m) for successive m from 1 to 16 
is made and /(m) is found to follow a power-law formula 
(inset of Fig. [2fb) ) . Having known the expression of f{m) , 
Eq.(2) can be used to determine Tc (Fig. [2fb)). The ma- 
jor error of this method results from the finite number 
of the fitted data points. The error can be evaluated in 
the following way. Consider that we have a sequence of 
data points of successive m up to m = m/, the fit for 
these data points gives Tc ^ . Then the error is given by 

iTr'-Tci= iTr^-rr^+Vrr^^^ -t^^^^^...^t^-Tc\ 

< J^ \T^ — T^^^\dx. We calculate T^^^ for successive 

rrif and find that |Tc^^ — TJ^^ | follows a power-law for- 
mula. In the case of Fig. [2fb), for example, the formula 
approximately equals O.SttiT^. With ruf = 20, the error 
is estimated to be smaller than 10~^. Actually our cal- 
culation for each Tdk) indicates that the error caused by 
the proposed finite size scaling method is no more than 
O(10~^). Detailed results of Tcik) are reported in the 
legend of Fig. [3l 

To find a general correlation of critical point and the 
network's structure, Tc is plotted as a function of the av- 
erage degree /c(Fig. [3]). It gives 



Tc{k) 



(3) 



where r = 0.49 ± 0.02. Rather interestingly, Tc(k) still 
scales as a power law but with a different exponent — r ?^ 
—0.5, in contrast to —1 found in ordinary site percolation. 
What is more nontrivial is that a general expression of 
tc{m^k) can be directly deduced from Eq.(3). Since for 
all 1 < m < A^, tc{m^k) is bounded by two power-law 
functions, i.e. Eq.(3) and tc{m = 1,/c) = 1/k (the critical 
point of ordinary percolation) , according to the argument 
in Ref. [26] critical point tc{m^ k) for all m must scale as a 
power law with k and its exponent depends on the specific 
value of m. 

tc{m,k)^k-^^'^\ (4) 

Fig. [3] demonstrates the validation of Eq.(4). The expo- 
nent X{m) is found to be well approximated by an arctan 
function 

A(m)^^ ^^TT- (5) 

^ ^ 4 arctan mO-4 ^ ^ 

Although the validation of Eq.(5) needs to be further ex- 
amined since it is a conjectured result, it coincides with 
A(l) = 1, X{m = N ^ (X)) = 0.5 and also fit the body well 
as shown in Fig. ^h). 

What differs from the ordinary percolation is that the 
scahng (Eq.(4)) is not always vahd but deviates from 
power law for small k when ?n > 1, as demonstrated in 
the inset of Fig. [^b) for the extreme case of m = N (i.e. 
Tc). Actually the deviation becomes more apparent as m 
grows. This deviation from the power law origins from two 
aspects: the unconnectedness of the underlying ER ran- 
dom graph for small k and the preferential selection for 
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Fig. 3: (a): Correlation of tc(m,k) with k for different m. All of 
them follow power-law correlation with exponent decreasing with 
m, ranging from r ~ 0.5 to 1. Blue line represents the result of 
ordinary percolation, i.e. tc(l,k) = 1/k. Brown line represents 
Eq.(3). Tc{k) (brown left-triangle) for /c = 4, 5, 6, 8, 10, 14 is 0.757(4), 
0.682(2), 0.633(3), 0.554(5), 0.484(3), 0.421(3), respectively, where 
the number in the bracket refers to the error of the last digit. Eg. 
0.757(4) means 0.757 ± 0.004, respectively. Note that this error is a 
total error, i.e. it includes errors from a finite-size scaling of tc(m) 
and Eq.(2) 

(b): Correlation of the exponent A(m) with m. Blue solid 
point is the simulation results while red line represents 
Eq.(5). Inset: Full correlation of Tc with k. The predic- 
tion by Eq.(6) (red line) coincides well with the simulation 
(blue solid point). In both simulations, numerical calcu- 
lations for tc{m^ /c, N) is carried out up to network size of 
10^ nodes and is averaged over 50 realizations. 



isolated nodes in explosive percolation. When /c = 2, for 
example, despite of the existence of the giant component 
of the underlying graph there are still about 20 percent 
of nodes forming isolated small clusters. In contrast to 
an absolutely uniform occupation in ordinary percolation, 
these isolated nodes are preferentially occupied before crit- 
icality in the explosive ones but have no contribution to 
the emergence of the giant percolated cluster. However, 
in a connected graph every node helps to bring out gi- 
ant cluster during the process of percolation. Therefore 
the constituent of tc{m > 1,/c) of small k includes many 
idle nodes and thus causes the deviation of power law. If 
m grows, the degree of preferential occupation of isolated 
nodes increases, so the deviation becomes more apparent. 
To make the explanation more convincing, Tc for small k 
will be reconstructed according to the above explanation. 
In the case of m = A/" ^ oo, it is expected that i) all the 
underlying isolated nodes will be occupied before critical- 
ity; ii) the whole network will not be percolated until the 
underlying giant component is percolated. Therefore Tc 
can be written as the sum of two parts: the relative size 
of those isolated nodes and the relative size of the occu- 
pied nodes at criticality within the giant component of the 
underlying graph, as expressed in the following equation 



Tc{k) = l-Sr^Tc{^{k))Sr 



(6) 



where Sr is the relative size of the giant component of the 
underlying graph which satisfies the well-known equation 

Sr + e-^^- = 1 [271. Tc{n{k)) is the critical point of the 
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percolation on the underlying giant component whose av- 
erage degree is n{k). Clearly 1 — 5^ is the relative size of 
the isolated nodes while Tc{n{k))Sr represents the contri- 
bution of the occupied nodes within the underlying giant 
component. Since ^r ^ as /c ^ 1, it is clear that the 
critical point occurs at 1 for /c = 1. For k > 1 the predic- 
tion of Tc depends on n{k) which is found numerically as 
n{k) = 2.337e-0-^^^^ + k. Therefore to obtain Tc{k), we 
need to know Tc{n{k))^ and consequently Tc(/^2(^)), where 
n2{k) = K.{K.{k))^ and so on. Finally, we will arrive at a 
certain / so that hii{k) > 4. Having known the validation 
of Eq.(3) for /c > 4, the above-mentioned method allows 
a recursive estimation for Tc. As shown in the inset of 
Fig. [3l^b), the prediction coincides the simulation results 
well, which confirms our explanation. 

Eq.(3) to Eq.(5), as the main contribution in our study, 
not only generalize some classical conclusions but could 
also give some new insights to the network robustness un- 
der attack. While most related studies designed their at- 
tack strategy based on some topology properties such as 
degree [2^ 125] , a practical attack strategy depends totally 
on how much the attacker knows about the network. If the 
attacker gets no information on network, random attack 
is the only option (the case of m = 1). But if the attacker 
knows the network structure every detail, he/she can find 
the best way (not just based on degree or some other prop- 
erties) to destroy the network (the case of m = N). Thus 
the practical meaning of m could be a measurement of the 
information that the attacker has on the network. Our 
study indicates that it is possible to understand the net- 
work robustness to any kind of attack based on any such 
information. However, since assessments of these informa- 
tion are usually impractical, Eq.(3) shows its great signifi- 
cance because Tc can be served as a warning that possible 
breakdown could ever occur once the attack process ex- 
ceeds this point. These concepts and ideas are nontrivial 
in the sense that they indicate a possible general frame- 
work for the study of network robustness. Indeed from the 
viewpoint of information, the network robustness might be 
defined as a strategy-independent function, namely tc{m) 
while previous strategy-dependent work can hardly pro- 
vide such general definition due to their unquantifiable 
model. 

Now let us turn to the discontinuity of the explosive 
site percolation. The explosive nature can be established 
in the following way [9]. Let to denotes the last step for 
which Si < A/'", and ti the first step for which Si > cN, 
where a and c are two positive constants. If there exists a 
group of a and c so that the interval At/N = {ti — to)/N 
vanishes in the thermodynamic limit, the percolation is 
said to be explosive. In the ordinary site percolation, such 
a and c never exist, leading to At always of order of TV. 
According to the above definition, percolation with zero 
critical point can never be explosive because the order pa- 
rameter cannot have macroscopic increase until sufficient 
(macroscopic) number of nodes are occupied. Since Eq.(3) 
indicates that for any dense ER network {k{N) -^ oo when 
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Fig. 4: (a): The scaling of At with A^. For m = 2, we find 
a ^ 0.5 and c ^ 0.15, leading to At ~ ArO-84 ^^^^^ j-^^gj p^^. 
m = 3, we find a ^ 0.6 and c ^ 0.25, leading to At -^ A^^-^^ (black 
line). Both of the interval At/N vanish in the thermodynamics limit. 
But for the traditional percolation, the interval always extends, i.e. 
At ~ 0(N), as indicated by green line, (b): The power-law relation- 
ship of AS'macc(the size of the largest jump of the order parameter) 
and network size N. The exponents for the case of m = 2, 5, 20 are 
measured as 0.15,0.088,0.05, respectively. All the results are ob- 
tained under network size up to 10^ nodes and are averaged over 50 
realizations. 



TV -^ oo) the critical point is zero regardless of tti, perco- 
lation on dense ER network must be non-explosive and 
continuous. Sparsity is thus a necessary condition for a 
percolation to show explosive nature. If k is finite, for 
m > 1 the percolation seems immediately becoming explo- 
sive with At ~ o{N) (Fig. Hl^a)). However this simulation 
result can be deceptive and misleading, as indicated by Ref 
J2T] . Indeed further studies on the largest jump of the or- 
der parameter indicates that the jump size decreases with 
TV as power law and vanishes when N -^ oo(Fig. Hl^b)). 
Thus the site explosive percolation is still continuous in 
the thermodynamic limit for these cases. 

To find a general conclusion for all m, finite size scaling 
analysis of S for various m is performed [2)fT2] . The theory 
of finite size scaling tells us that the order parameter obeys 
the relation S = N-^^''F[{t-tc)N^^''], where p and v are 
two critical exponents and F[*] is a universal function. 
If the percolation is discontinuous, this scaling relation 
trivially applies with /3 = 0. Otherwise studies of 5 as a 
function of the system size A^ yield the exponent of j3. The 



analysis finally gives the relation j3(m) 



m 



^(Fig.ETa)), 



indicating /3 ^ only when m ^ oo. The finite size 
scaling analysis is also applied to the susceptibility x which 
obeys x — N'^ ^ ^ G[{t — tc)N'^^ ^]., where 7 is another critical 
exponent and G[*] is the related universal function. The 
definition of x along with the scaling behavior of S and x 
give the relationship 7/1/ + P/v = 1, which can be used 
to check our measured value of the critical exponents [12] . 
As shown in Fig. [5fb), the sum {j3 + 7)/?^ is always 1 
with good approximation. Other critical exponents are 
also reported in Fig. [IJb). The power-law decay of j3{m) 
indicates that explosive site percolation on ER network 
is continuous for any finite m but could be discontinuous 
when m approaches infinite. 

In summary, the critical points of a class of explosive site 



p-4 



Criticality and Continuity of Explosive Site Percolation in Random Network 




Fig. 5: (a): Finite size scaling analysis of the explosive percola- 
tion for the case of ttt, = 3. Critical exponents can be determined 
by making all the data points of different network size collapse to- 
gether. For m = 2, 3, 5, 10, 20, (3 ^ 0.41, 0.26, 0.15, 0.074, 0.031 while 
^ = 1 for 777, = 1 as known for traditional site percolation. Inset: The 
scaling relationship of critical exponent (3 and m. (b): The relation- 
ship of the critical exponents and m. For m = 2,3, 5, 10, 20, (3/^ ~ 
0.221,0.164,0.105,0.055,0.025 and 7/1/ ^ 0.77,0.84,0.9,0.95,0.98. 
The sum ((3 + 7)/^ is always 1 with good approximation. All the 
results are obtained under network size up to 10^ nodes and are 
averaged over 50 realizations. 



percolation in ER random network are found to scale with 
the average degree and their exponents range from —1 to 
— r ^ —0.5. Heuristic discussions are made to uncover 
their possible implication to network robustness. Note 
that although these results are obtained under our spe- 
cific percolation model, they could be even more general. 
Since the explosive percolation, as the model proposed 
here, controls directly the growth of the giant cluster, it is 
expected that when m = N^ the critical point is delayed 
to the greatest extent compared to any other possible site 
percolation (i.e. the upper bound provided by Eq. (3) may 
be valid for any site percolation on ER network). If this 
assumption is true, for any site percolation satisfying crit- 
ical point tc > l/k on ER network, tc must follow Eq. (4) 
because such critical point is bounded by l/k and Eq. (3). 
Thus Eq. (4) could be a general and model-independent 
result. The explosive nature, according to its definition, 
is deduced to occur only on sparse network. So sparsity is 
a necessary condition for the explosive site percolation to 
be discontinuous. When this condition is fulfilled, analy- 
sis on the critical exponents by finite size scaling method 
indicates a power-law decay of j3{m). Other critical expo- 
nents are also presented to show their consistence. These 
results reveal that discontinuity of the class of explosive 
site percolation could happen only in sparse networks and 
only when m ^ 00, while for any finite m the explosive 
percolation is still continuous. 
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